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performed. At tree level, only charged pions couple to the initial state photons and 



CN ' A quark model calculation of the processes 77 — > tt + tt and 77 — ► 7r 7r° is 
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00 ' neutral pions are not expected in the final state. However, a small but significant 
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77 — > 7r°7T cross section is observed. We demonstrate that this may be accounted for 
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by a rotation in isospin space induced by final state interactions. The resulting 7r + 7r~ 

Oh! 

O 1 cross section is in good agreement with experiment while the 7r°7r° cross section is in 

cu : 

■ qualitative agreement with the data. 

X 

.«?. 

I. INTRODUCTION 

The constituent quark model has proven exceptionally successful in describing the prop- 
erties of hadrons 0. Nevertheless it is universally acknowledged that the quark model 
spectroscopic analysis is overly simplistic, in part because couplings of resonances to other 
hadronic states are not included. In particular, the quark model at its most basic level does 
not include final state interactions and the effects of coupling resonances to decay chan- 
nels |U . One can view the constituent quark model as representing the lowest components of 
a Fock space expansion while final state interactions and coupled channel effects represent 



higher order contributions. It is possible for these contributions to shift significantly the 
resonance pole position. 

In this paper we study the contributions of final state interactions to the process 77 — ► tttt 
for relatively low center of mass energies. In this energy regime the photons have sufficiently 
long wavelength that they do not resolve the pions into their constituent quarks. This asser- 
tion has been formalized in the Low theorem ||, which states that scalar QED is adequate 
to describe 77 — > 7r + 7r~ near threshold. However, scalar QED cannot explain the observed 
finite cross section to neutral pions. We will find that including final state interactions gives 
qualitative agreement with the data. This is nontrivial in view of the subtlety of the re- 
action, the expected importance of relativistic effects, and the simplifying assumptions we 
have made. Essentially, the 7r°7r° final state arises from the difference in the 1=0 and 1=2 tttt 
potentials which contribute to the rescattering of the final state pions in 77 — > ttti. Although 
the final state interaction contributions have been studied previously in the quark model 
||, that initial study used the Fermi approximation to final state interactions, did not take 
into account the differences in the 1=0 and 1=2 potentials (which are crucial to describing 
the 7r°7r° final state), and employed rescaled pion potentials. As a result, the model of Ref. 
U predicts a substantial enhancement in charged pion production and fails when applied to 
neutral pion production 0. Our conclusions differ because we have avoided these pitfalls. 

Before proceeding we note that this process has been studied in other approaches which 
give good agreement with experiment. In particular, a recent calculation using chiral per- 
turbation theory to two loops || agrees with the data as do calculations using dispersion 
relations with phase shift data from tttt scattering and constraints from unitarity, analytic- 
ity, and crossing 0. Our purpose here is therefore not so much to explain the data but to 
expand the regime of applicability of the quark model and to gain some experience so that 
we may apply these methods to situations where the techniques just mentioned are inappli- 
cable. Indeed, a large body of data exists on vector-vector production from 77. Much of 
this data is poorly understood and has led to wide spread speculation |TD[ . We believe that 



the results presented here will aid in clarifying this enigmatic situation. 



We begin in section II by summarizing the scalar QED results for 77 — > 7r + 7r~. In 
particular, we expand the cross section in terms of partial wave amplitudes. The expressions 
for the cross sections including final state interactions are derived in section III. In section 
IV we present our results and, finally, a summary is given in section V. 



II. THE SCALAR QED CROSS SECTION FOR 77 -> vr+vr 

At low energies, long wavelength photons do not probe the quark constituents of pions. 
Thus the pions may be rigorously treated as point particles interacting via scalar QED || . 
The cross section for the process 77 — > 7r + 7r _ may be derived from the amplitudes shown in 
Fig. 1 and is given by 
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where x = 4m 2 fs, s = E 2 m is the Mandelstam variable, and m is the mass of the pion. 
Partial waves amplitudes, which are required when considering final state interactions, are 
obtained from the helicity amplitudes: 
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where 8 and are the spherical polar coordinates of the outgoing 7r + and the subscripts 
+ and — are the helicities of the two incoming photons. Partial wave amplitudes may be 
extracted as follows: 
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One finds that f£jfc L = f LMh is only non-zero for M L = and that f£ ML = f L \. lL is only 
non-zero for Ml = 2. The cross section is then written in terms of the partial waves as 
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where S = 1/(1 + S nn ) is a statistical factor required for the 7r°7r° channel (to be discussed 
below). The necessary partial waves are given by: 
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III. FINAL STATE INTERACTIONS 

Final state interactions (FSI) are additional interactions between particles in the final 
state that are not included in the basic interaction. In our problem the basic interaction is 
the scalar QED process described in the previous section and the FSI's are purely hadronic 
and are described by effective potentials between the outgoing particles. We use the distorted 
wave Born approximation in our calculations 0, which relies on nonrelativistic quantum 
mechanics. We include relativistic phase space where possible and use the relativistic scalar 
QED result for the basic interaction. This mixture of nonrelativistic and relativistic elements 
is typical of the quark model. 

The first step in including FSI's is to find the wavefunction describing the relative motion 
of the final-state particles due to the FSI potentials alone. We choose to do this by solving 
a two body Schrodinger equation with an effective irir interaction. The scattered waves are 
written in terms of spherical waves as 
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where ui(k,r), the radial solution of the Schrodinger equation, has the asymptotic form 

ui(k,r) r - z ^' -j- sin (Ax - ±/7r + (7) 

Here /c is the relative momentum between the particles, k 2 = 2fiE, where \x is the reduced 
mass of the particles and E is the total energy of the system. 



The 7T7T potentials used to derive the final state wavefunction are taken from Refs. []TT1 , [T2 
These are obtained by calculating the perturbative T-matrix for a given process in the 
nonrelativistic quark model. This T-matrix was then equated to a model pointlike T-matrix 
to obtain an effective interaction. The quark model calculation includes one gluon exchange 
and confinement interactions followed by quark rearrangement into final state color singlets. 
For the 1=0 channel there is an additional contribution from qq annihilation. The resulting 
effective potentials are parametrized in the form 

V^(r) = V e' h (^) . (8) 

The parameters, Vo and tq, are reproduced in Table 1. FSI corrections were only made to 
the S-wave amplitudes because corrections to D-waves were found to have a negligible effect 
on the cross section. 

Final state interactions are incorporated with the standard two-potential relation [Q,[3|]: 
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where ip^*(k) is the complex conjugate of the momentum-space solution of the Schrodinger 
equation with the FSI potential V nn , W is the basic electromagnetic interaction, and 0£ is 
a plane wave. The corrected QED amplitude is 
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where M(k) is the uncorrected scalar QED amplitude, ji(kr) is a spherical Bessel func- 
tion and the yfs factors are introduced to relate the relativistic QED amplitudes to the 
nonrelativistic amplitudes used in the calculation. 

In order to include the FSI correction for each partial wave we decompose the amplitudes 
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where the corrected partial wave amplitudes are given by 
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Helicity labels have been omitted. In the absence of an FSI potential, 5l — > and 
u L (kf,r) — > j L (kfr), the orthogonality of the spherical Bessel functions yields a delta func- 
tion, and we recover fi,M L (s(kf)) — > fLM L (s(kf)) as expected. 

The relationship between the charged pion basis of the physical scalar QED reaction and 
the isospin basis of the potentials is 

k + = JI |00> + Ji|20>, 



k = -y-|oo> + y-|2o>. 

Thus, the physical pion wavefunctions can be written in terms of the momentum space 
solutions of the Schrodinger equation in the isospin basis: 
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The corrected partial waves for it + it production are thus given by 
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with an analogous expression for the 7r°7r° final state. Here S ! L is the nn phase shift in the 
isospin /, partial wave L channel. To simplify the following expressions we define the isospin 
amplitudes 
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This yields the final state enhanced probability for 7r + 7r production: 
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and for 7r°7r° production: 
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With no nn potentials (7la/ l —* Ilm l and 5£ = and the original cross sections are recovered. 

The last step is to account for the limited polar acceptance of the nn experimental data. 
Integrating over a finite solid angle (— cos# acc < cos# < cos# acc ) to obtain the uncorrected 
cross section for limited polar acceptance is straightforward: Eq. [I] is replaced by 
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However, the expression in terms of partial waves is affected more drastically as the 
spherical harmonics are no longer orthogonal to each other and the contributions of the 
partial waves to the total cross section can no longer be separated. The expression replacing 
Eq. |] becomes an infinite sum. 
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Because the L = partial wave is the only one affected significantly by FSI's, we assume 
that all of the observed 77 — > 7r°7r° events are in an L = state. Since that distribution is 



spherically symmetric, we can correct the data by simply dividing by the value of cos# acc 
appropriate to the experiment in question. 

For 77 — > 7r + 7r~ the total observed cross section is a mixture of all partial waves. We 
again assume that only the L = wave is affected by FSI's. To simplify our calculations we 
define the harmonic integrals 
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and note that h™ l2 = and hj ll2 = hf ll2 (h is independent of m).. We can then write 
the total cross section with only the L = partial wave corrected for FSI's as the total 
uncorrected cross section, plus an infinite number of correction terms, each of which is the 
difference between the corrected and uncorrected values of a term in Eq. |17| Only those 
terms involving / 00 , which goes to / E s 7 |: +7r _, need correcting. Eq. 4 is replaced by 
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Excellent convergence over the relevant energy range is obtained when the series is truncated 
after L = 6. Fortunately, the three experiments whose data we use for comparisons all have 
the same limited polar acceptance (| cos^l < 0.6), so the results need be corrected for only 
one value of cos6L rf .. 



IV. RESULTS 

Equation (13) demonstrates that the isospin dependence of the nn production amplitudes 
is generated through the isospin dependence of the final state interactions. The fact that 
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the isoscalar and isotensor tttt interactions differ markedly causes a rotation in isospin space 
which is reflected in a nonzero 7r 7r° amplitude (see Equation (15)). Indeed, strong, but 
isospin independent, nir interactions would not give rise to neutral pions in the final state. 

Equation (15) reveals that the details of the neutral pion prediction depend on the 
scalar QED production mechanism, the isoscalar and isotensor hadronic phase shifts, and 
FSI-induced modifications to the outgoing pion-pion wavefunction. We shall assume that the 
tree order scalar QED amplitude captures most of the production mechanism (we note that 
direct 77 coupling to the fo and rescattering to 7r°7r° is possible - this is discussed below) 
and that the joining of relativistic and nonrelativistic formalisms has not overly distorted 
the predictions. 

The S-wave elastic phase shifts produced by the irir potentials are shown as solid lines in 
Figures 3 (/ = 2) and 4 (/ = 0). It is clear that both underestimate the strength of the nir 



interaction. However, it has been shown [17| that the majority of this discrepancy may be 



alleviated by employing relativistic kinematics in the description of the tth phase shift. This 
yields the dashed curves in Figures 3 and 4. The agreement in the isotensor sector is quite 
good. The isoscalar prediction is still somewhat low; however, the isoscalar sector differs 
from the isotensor sector in that qq annihilation will lead to resonance coupling. Indeed, 



the authors of Ref. |Tl| have shown that very good agreement may be achieved once the 
neglected effects of mixing with the / (980) and / (1350) are included. Since resonance 
coupling will also induce wavefunction effects in our problem, we have chosen to examine 
the importance of the phase shifts by simply fitting them. The fits are indicated as dotted 
lines in Figures 3 and 4. The differences between the predictions for neutral pion production 
for the relativized model (dashed lines) and the fits (dotted lines) will be examined below. 

The final contribution to the 7r°7r° cross section came from the FSI-induced distortions 
to the pion-pion wavefunction. We shall assume that this is accurately described by our 
potential formalism. This may appear doubtful because the potential phase shifts did not 
agree well with elastic 7T7T scattering data. However, as noted above, the majority of the 
discrepancy is due to the nonrelativistic phase space necessitated by the numerical solution 



of the Schrodinger equation, rather than the potential itself. We are thus confident that the 
effective potentials can describe the wavefunction distortion. However, the model neglects 
additional distortions to the wavefunction caused by transitions to virtual resonance states 
(the /o). For the moment, we will assume that these are small near threshold since the 
/o(980) is narrow and the /o(1350) is distant. Indeed, calculations with a relativistic Breit- 
Wigner yield the following cross sections for 7r°7r° at threshold: 0.015 nb for /o(980), 0.00005 
nb for / (1275), and 0.15 nb for / 2 (1275). 

Figure 5b shows the predicted cross section for 77 — > 7r + 7r - along with data from three 
experiments. There is very little difference between predictions, indicating the lack of im- 
portance of FSI effects in the charged channel. The Born result agrees fairly well with the 
data. The Born result also agrees well with the dispersion analysis of Ref. |J for y/s < 350 
MeV, where there is no data. 

The dashed line of Figure 5a is the predicted cross section for 7r°7r° production which 
is based on the relativized tttc phase shift predictions of Ref. ]I2] and a nonrelativistic titt 
final state wavefunction. Recall that this curve therefore does not incorporate resonance 
effects and does not fit the isoscalar tttt scattering data well. The result presented in Figure 
5a agrees well with data near threshold - an indication that the near threshold neutral 
pion data can be explained by the postulated mechanism of final state interaction-induced 
isospin rotation. Nevertheless, the dashed curve falls below the data by a factor of three or 
four above energies of 350 MeV. It is therefore of interest to test the effects of intermediate 
resonance states. We model this by fitting the experimental S-wave tctt phase shifts as 
described above. Note that we are ignoring wavefunction distortion caused by dynamical 
resonance effects. The resulting prediction is given by the dotted line in Figure 5a. The 
agreement with data is substantially improved. This is due to the increased accuracy of the 
description of the 1 = scattering data and is a strong indication that resonance effects (as 
manifested in the 1 = phase shift) can be important even near threshold. 

The agreement with experiment is not perfect - we believe that this is likely due to 
explicit dynamical resonance effects in the final state scattering wavefunction. It should be 

10 



possible to test this assertion by including resonances in a coupled channel description of 
the final state interactions. Nevertheless, the main point is established - it is possible to 
describe the rather complicated processes leading to a 7r°7r° final state with a sufficiently 
detailed analysis of the production mechanism and the final state interactions. 



V. CONCLUSIONS 

The reaction 77 — ► 7r°7r° presents an interesting theoretical challenge. Neutral pions 
do not couple directly to low energy photons so that direct techniques based on the Low 
theorem cannot be used to describe the cross section. Calculations tend to be difficult; 
this is illustrated by efforts in chiral perturbation theory which have found it necessary to 
introduce three new counterterms, work to two loops, and to unitarize. We have presented a 
relatively simple explanation of the data. In our model neutral pion production proceeds via 
charged pion production followed by final state hadronic interactions. The it + tt~ initial state 
is described well by the Born order scalar QED amplitude, in accord with Low's theorem. 
This amplitude is then convoluted with the final state pion-pion wavefunction to arrive at 
a prediction for 7r°7r° production. 

The result for 77 — > 7r + 7r~ is very insensitive to final state interactions - a result that 



follows from isospin algebra and the weakness of the final state interactions p8 |. However, 
the prediction for 7r 7r depends crucially on the interference between 1 = and 1 = 2 final 
state interactions, and is therefore driven by the well known differences in these processes 
-1 = 2 S-wave tttt scattering is repulsive while I = is attractive. It was shown in Ref. 
fll| that roughly half of the scattering strength at low energy in / = is due to the / (980) 
and /o(1350) resonances - even at threshold. This leads to the surprising conclusion that 
resonance coupling is important to getting the magnitude of the 7r°7r° cross section correct 
near threshold. 

This calculation was undertaken with an eye to extensions to heavier systems. In partic- 
ular, 77 production of vector-vector states has a long, and enigmatic, history. For example, 
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the large ratio (7(77 — > p p°)/°"(77 — > p + p~) is problematic and has engendered much 
speculation. This is due, is part, to the lack of clear theoretical guidance - chiral per- 
turbation theory and dispersion methods are inapplicable for these processes. Of course 
such constraints do not exist on the quark model and we look forward to applying it to an 
explication of this interesting area of hadronic physics. 
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FIGURES 

FIG. 1. The tree-level Feynman diagrams for the scalar QED contributions to 77 — ► it + it~ . 
The labels denote the particle momenta (pi) and photon polarizations (e$). 

FIG. 2. The ttit potentials vs r for / = 0, L = (solid line), I = 0, L = 2 (dashed line), 1 = 2, 
L = (dotted line), and I = 2, L = 2 (dot-dot-dashed line). The potentials are given by Eqn. || 
with the parameters given in Table I. 

FIG. 3. The 1 = 2 tttt scattering phase shift vs. y/s using ttit potentials (solid line), the 
quark model prediction (dashed line), and a linear approximation to the experimental data, 
djzl = -0.00062 rad/MeV x (y/s- 250 MeV) (dotted line). The data are from Hoogland et al. 
(circles) |]l3| - we show the results of both of their methods; and Prukop et al. (triangles) |l4j 

- we show the results of their first fit. Some of the data shown in this figure were obtained from 
References flj| and ]r|. The horizontal error bars on the data indicate bin size; for the vertical 
error bars all of the given errors were added in quadrature. 

FIG. 4. The 1 = ttit scattering phase shift vs. \/s, using tttt potentials (solid line), the 
quark model prediction without resonance contributions (see text) (dashed line), and a linear 
approximation to the experimental data, <5{z.q = 0.0027 rad/MeV x (y/s — 250 MeV) (dotted line). 
The data are from Mukhin et al. (circles) Rosselet et al. (inverted triangles) [[H]] - the 

horizontal bars only approximate their bins, and their data is actually for <5§ — 5\ - we have used 
the 8q data extracted from it by Li et al. Estabrooks and Martin |^0| - we show the results 
of both their s- (squares) and t-channel fits (diamonds); and Protopopescu et al. (triangles) |2l| 

- we show the results of their case 1. For further comments see Fig. 3. 
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FIG. 5. (a) (7(77 — > 7r°7r°) vs y/s, using phase shifts from the quark model (dashed line), from 
linear fits to the data (dotted line), and experimental data. The data are from Bienlein (et al. ) 
(closed circles) p^| , Marsiske et al. (open circles) [£3| and Edwards et al. (closed triangles) f24fl . 
The data have been corrected to full polar acceptance, |cos#| < 1.0. (b) 17(77 — ¥ tt + tt~) vs. 
\/s with the same line labelling as (a). The curves have been corrected to have a limited polar 
acceptance to match the data: | cos#| < 0.6. The data are from Behrend et al. (closed circles) p5| , 
Boyer et al. (open circles) ]26| , and Aihara et al. (closed triangles) |27j. For additional comments 
see Figure 3. 



16 



Table 1: The parameters of the 7r — -k potentials used in this work [17,11]]. The t-channel gluon 
exchange potentials have two contributions, due to colour-hyperfine and confinement terms, which 
must be summed. The potentials use the parametrization of Eqn. with Vq in GeV, tq in GeV -1 . 



Pion State 



L 






2 
2 




2 

2 



t-channel gluon exchange 



hyper fine 



Vn 



-0.392 
0.786 

-0.044 
0.088 



r 



1.36 
1.36 
1.40 
1.40 



s-channel 
confinement Igluon exchange 



Vn 



-0.024 
0.047 
0.0175 
-0.035 



r 



2.29 
2.29 
1.49 
1.49 



-1.75 



r 



1.48 
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fa) rt>) (c) 

Fig 1: The tree-level Feynman diagrams for the scalar QED contributions to 77 — ► tt + it~. The 
labels denote the particle momenta (pi) and photon polarizations (e^). 




-2500 



r (fm) 

Fig 2: The -kit potentials vs r for J = 0, L — (solid line), I = 0, L = 2 (dashed line), 
/ = 2, L = (dotted line), and I = 2, L = 2 (dot-dot-dashed line). The potentials are given 
by Eqn. |8] with the parameters given in Table I. 
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Figure 3: The 1 = 2 tttt scattering phase shift vs. yfs using tttt potentials (solid line), the 
quark model prediction [O] (dashed line), and a linear approximation to the experimental data, 
djzl = -0.00062 rad/MeV x (y/s - 250 MeV) (dotted line). The data are from Hoogland et al. 
(circles) |]l3| - we show the results of both of their methods; and Prukop et al. (triangles) |l4j] 
- we show the results of their first fit. Some of the data shown in this figure were obtained from 
References [15] and ]l(|. The horizontal error bars on the data indicate bin size; for the vertical 
error bars all of the given errors were added in quadrature. 
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Figure 4: The 1 = ir-ir scattering phase shift vs. y/s, using -kit potentials (solid line), the 
quark model prediction without resonance contributions (see text) JT1J (dashed line), and a linear 
approximation to the experimental data, SjzPq = 0.0027 rad/MeV x (y/s — 250 MeV) (dotted line). 
The data are from Mukhin et al. (circles) |l8|| ; Rosselet et al. (inverted triangles) jj/J] - the 
horizontal bars only approximate their bins, and their data is actually for 5q — 5\ - we have used 
the 5q data extracted from it by Li et al. fll| ; Estabrooks and Martin pfj - we show the results 
of both their s- (squares) and t-channel fits (diamonds); and Protopopescu et al. (triangles) pi] 
- we show the results of their case 1. For further comments see Fig. 3. 
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Figure 5: (a) 17(77 ~^ vr°7r°) vs ^/s, using phase shifts from the quark model (dashed line), from 
linear fits to the data (dotted line), and experimental data. The data are from Bienlein (et al. 
) (closed circles) |2^| , Marsiske et al. (open circles) J23|, and Edwards et al. (closed triangles) 
pif . The data have been corrected to full polar acceptance, | cos#| < 1.0. (b) 17(77 — * vr + 7r _ ) vs. 
1/s with the same line labelling as (a). The curves have been corrected to have a limited polar 
acceptance to match the data: | cos#| < 0.6. The data are from Behrend et al. (closed circles) [25] , 
Boyer et al. (open circles) f26fl , and Aihara et al. (closed triangles) J27J. For additional comments 
see Figure 3. 
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